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for  Prediction  of  UV  Radiation  from 
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Numerical  study  of  a  5,900  lb  solid  propellant  thruster  plume  at  an  altitude  of  118  km 
is  performed  using  a  combined  multi-step  continuum/kinetic  approach.  The  Navier-Stokes 
equations  are  solved  for  the  flow  inside  the  nozzle  and  the  first  several  meters  into  the 
plume.  The  DSMC  method  is  used  to  simulate  the  remaining  plume  and  the  plume-free 
stream  interaction.  A  Monte  Carlo  based  radiation  code  is  applied  in  an  overlay  mode  to 
calculate  UV  radiation  in  the  near  and  far  field  of  the  plume.  The  computations  take  into 
account  both  alumina  particles  and  soot.  The  effect  of  alumina  particle  emissivities  and 
size  distribution,  as  well  as  soot  concentration,  on  UV  radiation  is  clarified.  Comparison 
of  numerical  results  with  available  UV  measurements  is  conducted. 


I.  Introduction 

The  success  of  missile  signature  analysis  depends  heavily  on  the  reliability  of  signature  predictions  for 
all  stages  of  the  boost  trajectory.  Prediction  reliability,  in  turn,  is  significantly  hampered  by  a  large  number 
of  uncertainties  and  unknowns  currently  associated  with  the  physics  and  gas  dynamics  of  rocket  propulsion 
system  exhaust  plumes  interacting  with  the  surrounding  atmosphere.  An  essential  feature  of  solid  propellant 
rocket  plumes  that  complicates  flow  analysis  is  related  to  the  formation  of  particulates  of  varying  types,  sizes, 
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and  loadings.  This  makes  rocket  plume  flows  two-phase,  with  one  phase  consisting  of  gas  atoms  and  molecules 
and  the  other  composed  of  nano-  and  micro-scale  particulates.  Disregard  of  the  two-phase  character  of  rocket 
plumes  usually  results  in  unacceptable  loss  in  prediction  accuracy.  The  specifics  of  interaction  between  the 
two  phases  and  its  impact  on  signatures  and  contamination  depends  on  the  propellant  composition  and 
thruster  operation.  The  main  particulate  types  are  soot,  propellant  droplets,  and  aluminum  oxide  particles. 
The  large  mass  fractions  of  aluminum  oxide  particles  have  a  significant  impact  on  the  gas  flow  inside  the  nozzle 
(and  therefore  thruster  performance),  the  near  field  plume  structure,  and  the  far  Held  plume-atmosphere 
interaction  structure.  The  strong  influence  from  these  particles  on  the  flow  structure  inside  the  nozzle  and 
in  the  plume  results  in  a  dramatic  effect  on  radiation  signatures. 

There  are  a  number  of  physical  processes  and  quantities  that  affect  radiation  from  rocket  plumes  in  the 
UV  region,  most  are  associated  with  various  particulate  properties  and  propellant /motor  parameters.  In 
the  related  previous  work,1  one  of  these  processes  was  considered  in  detail,  namely  photon  scattering  on 
alumina  particles  in  a  Star-27  motor  plume  at  118  km.  The  numerical  results  were  compared  with  available 
on-board  flight  measurements2  for  two  lines  of  sight;  near  field  (4°  from  the  plume  axis)  and  far  field  (25°). 
It  was  shown  that  the  main  process  determining  the  far  field  radiation  is  photon  scattering  on  submicron 
particles.  Accounting  for  this  effect  increases  the  predicted  radiance  by  a  factor  of  300,  and  the  value  of  the 
full  radiance  at  230  nm  in  this  case  was  approximately  a  factor  of  four  lower  than  the  measurements.  The 
predicted  near  field  radiation  was  only  about  30%  higher  than  the  measurements  at  230  nm.  This  difference 
increases  to  an  order  of  magnitude  at  400  nm. 

The  main  goal  of  this  work,  that  is  a  continuation  of  Ref.  1,  is  to  analyze  possible  processes  for  the 
difference  between  computed  and  experimental  results,  at  both  inband  and  spectral  levels.  Modeling  of 
these  processes  uses  engineering  judgement  to  define  reasonable  bounds  for  the  parameters  of  each  process. 
Similar  to  Ref.  1,  a  multi-step  continuum- kinetic  approach  is  used.  Continuum  Navier-Stokes  solvers,  VIPER 
and  CFD-I — b,  are  used  to  predict  the  flow  inside  of  the  nozzle  and  in  the  near  field  of  the  plume,  respectively. 
A  kinetic  DSMC-based  solver,  SMILE,  is  used  in  the  mid-  and  far  fields  of  the  flow.  The  final  flow  solution  is 
obtained  through  successive  application  of  these  solvers  from  the  inside  of  the  nozzle  to  the  outside,  with  the 
information  transferred  through  boundary  conditions  (starting  surfaces).  Finally,  a  Monte  Carlo  radiation 
code  NEMO  is  applied  to  calculate  UV  radiation. 
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The  following  important  factors  are  examined. 

(i)  The  impact  of  soot  particles  on  spectra  in  the  far  field.  Soot  particles,  common  to  both 
hydrocarbon  and  solid  propellant  systems,  do  not  significantly  affect  the  plume  gas  flow  due  to  their  small 
concentrations,  but  are  of  interest  due  to  their  possible  influence  on  radiation  signatures.  Soot  particles  are 
included  in  the  present  simulations  in  addition  to  the  alumina  particles  considered  previously,  and  the  effect 
of  their  loading  and  size  distribution  is  clarified  below.  Another  possible  source  of  small  particulates  in  the 
plume,  related  to  the  homogeneous  and  heterogeneous  condensation  processes,  is  not  considered  in  this  work. 

(ii)  Alumina  particle  emissivities.  There  are  a  number  of  particle  material  properties  that  may  po¬ 
tentially  impact  the  gas  and  particle  parameters  in  the  plume.  One  of  the  most  important  properties  related 
to  radiation  signatures  is  particle  emissivity.  Particle  emissivity  can  be  accurately  calculated  for  spherical 
particles  using  Mie  theory,  provided  the  complex  index  of  refraction  is  defined  as  a  function  of  temperature 
and  wavelength.  Comparison  of  the  complex  index  of  refraction  deduced  from  available  experimental  data 
shows  differences  of  up  to  two  orders  of  magnitude  observed  between  similar  experiments.  Since  the  real 
part  of  the  refraction  index  is  nearly  constant  and  much  larger  than  the  imaginary  part,  this  two  orders 
of  magnitude  difference  in  the  imaginary  part  translates  to  a  two  orders  of  magnitude  difference  in  particle 
emissivities.  In  the  present  work,  existing  theoretical  and  empirical  correlations  for  the  particle  refraction 
index  are  analyzed,  and  the  sensitivity  of  radiation  signatures  from  a  Star  27  motor  to  particle  emissivities 
are  examined. 

(iii)  Alumina  particle  size  distribution.  Experimental  evidence4  indicates  the  presence  of  a  sub¬ 
stantial  number  of  submicron  particles  in  solid  propellant  rocket  plumes.  The  size  distribution  of  submicron 
particles  exhibits  substantial  variability  for  different  motor  designs.3,4  A  submicron  particle  size  distribution 
for  the  Star-27  motor  has  not  been  measured  to  the  authors’  knowledge,  so  three  models  were  developed 
to  test  the  sensitivity  of  UV  emissions  to  the  assumed  size  distribution.  The  three  distributions  span  four 
orders  of  magnitude  in  assumed  mass  fraction  for  0.2  micron  diameter  particles.  This  range  of  particle  sizes 
tests  the  importance  of  small  alumina  particle  scattering  on  the  25°  model  data  comparisons. 
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II.  Flow  Conditions 


Availability  of  detailed  in-flight  measurements  of  UV  radiation  from  second  and  third  stage  plumes  of  the 
Strypi  XI  rocket  provides  a  solid  basis  for  numerical  model  testing  and  validation.5  These  measurements  were 
conducted  during  the  Bow  Shock  Ultraviolet  2  Experiment  flown  in  early  1991.  As  part  of  this  experiment, 
inband  and  spectrally  resolved  radiation  from  the  Antares  II  motor  (second  stage)  and  the  Star  27  motor 
(third  stage)  in  the  range  from  200  to  400  nm  were  recorded  for  flight  altitudes  from  approximately  104  km  to 
119  km.  For  data  collection,  aft-viewing  photometers  and  spectrometers  installed  on  extractable  periscopes 
were  used.  Two  lines  of  sight  were  examined,  4  deg  and  25  deg  from  the  rocket  centerline. 

In  this  work,  the  point  of  the  flight  trajectory  that  has  been  numerically  examined  corresponds  to  the  Star 
27  steady  state  operation  at  an  altitude  of  118  km.  The  geometric  setup  used  to  approximate  experimental 
conditions  is  shown  in  Fig.  1.  The  free  stream  temperature  and  density  are  330  K  and  5.21x  1017  molecule/m3, 
and  the  atmosphere  species  mole  fractions  are  66.7%  N2,  7.6%  02,  and  25.7%  O.  The  rocket  moves  with 
velocity  of  2850  m/s.  Zero  angle  of  attack  was  assumed  in  the  computations.  The  nozzle  thrust  is  about 
5,900  lb.  The  exhaust  gas  composition  is  listed  in  Table  1. 

Three  alumina  particle  size  distributions  are  considered,  all  of  them  combining  to  the  total  mass  fraction 
of  the  alumina  particles  in  the  exhaust  of  29.5%.  The  first  distribution  is  a  5  particle  bin  distribution  used 
in  the  previous  work,1  with  the  mass  fraction  of  subnricron,  0.2/mr  diameter,  particles  of  0.15  relative  to  the 
total  alumina  mass  fraction.  The  second,  14  particle  bin  distribution,  uses  the  same  combined  submicron 
particle  mass  fraction  of  0.15,  but  the  number  of  sub-micron  particle  bins  is  11  and  not  1.  The  third,  15 
particle  bin,  distribution  is  based  on  the  Worster6  particle  size  model.  The  Worster  distribution  produced 
the  smallest  mass  fractions  for  the  submicron  particles.  The  alumina  particle  mass  fractions  for  these  three 
distributions  are  shown  in  Fig.  2. 

Two  different  mass  fractions  of  soot  particles  were  considered,  0.3  and  3%  of  the  total  propellant  mass. 
The  distribution  of  soot  particle  sizes  was  represented  by  three  bins,  40,  100,  and  200  nm.  All  bins  are 
assumed  to  have  |  of  the  total  soot  mass  fraction.  This  soot  size  distribution  was  used  for  all  three  alumina 
particle  size  distribution  cases. 
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III.  Numerical  Approaches 


A  multi-step  numerical  approach  is  used  for  modeling  the  Star-27  rocket  plume.  Detailed  description  of 
the  approach  and  accuracy-related  issues  is  given  in  Ref.1  Only  a  brief  outline  is  presented  below. 

Step  1:  flow  inside  the  nozzle.  An  axisymmetric  Parabolized  Navier-Stokes  (PNS)  code  VIPER8,9 
is  used  that  includes  finite  rate  gas  chemistry,  multiphase  capability  (via  a  two-way  coupled  Lagrangian 
method),  and  a  variety  of  mostly  empirical  models  for  gas-particulate  interaction  and  particulate  evolution 
phenomena.  The  PNS  scheme  is  applied  from  the  sonic  line  near  the  throat  of  the  nozzle  to  the  exit  plane 
of  the  nozzle.  The  combustion  chamber  pressure  and  temperature  were  assumed  to  be  34.3  atm  and  3450 
K,  respectively. 

Step  2:  plume  near  field.  An  unstructured,  chemically  reacting,  multi-phase,  turbulent,  Reynolds 
Averaged  Navier-Stokes  (RANS)  code  CFD++10’ 11  is  applied  with  the  macroparameters  at  the  nozzle  exit 
obtained  at  Step  1  used  as  the  plume  inflow  boundary  conditions.  A  fully  implicit,  2nd  order  in  space, 
Harten,  Lax,  van  Leer,  Contact  discontinuity  (HLLC)  Riemann  approximation  algorithm,  and  a  modified  two 
equation  k-epsilon  turbulence  model  are  used.  An  Eulerian  based  Henderson  drag  coefficient  model12  is  used 
for  the  gas-particle  flow  interactions.  All  the  2D  axisymmetric  grids  used  a  baseline  mesh  of  approximately 
50,000  quadrilateral  nodes. 

Step  3:  plume  far  field.  A  DSMC  based  computational  tool  SMILE13  is  chosen  to  model  the  plume  far 
held,  with  a  starting  surface  at  about  2  m  from  the  nozzle  exit  plane  obtained  using  CFD+- 1-.  The  majorant 
frequency  scheme14  was  used  to  calculate  intermolecular  interactions.  The  intermolecular  potential  was 
assumed  to  be  a  variable  hard  sphere.15  Energy  redistribution  between  the  internal  and  translational  modes 
was  performed  in  accordance  with  the  Larsen-Borgnakke  model.  Temperature-dependent  relaxation  numbers 
were  used.  To  achieve  adequate  numerical  accuracy,  multi-domain  DSMC  computations  were  performed,  with 
combined  44  million  simulated  molecules  and  12.5  million  collision  cells. 

To  analyze  the  impact  of  the  alumina  particle  emissivity  and  size  distribution,  four  different  cases  are 
considered  for  the  flow  outside  the  nozzle, 

Case  1.  5  bin  distribution  of  alumina  particles  with  emissivities  computed  using  Ref.  16, 

Case  2.  5  bin  distribution  of  alumina  particles  with  emissivities  computed  using  Ref.  17, 

Case  3.  14  bin  distribution  of  alumina  particles  with  emissivities  computed  using  Ref.  17, 
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Case  4.  15  bin  distribution  of  alumina  particles  with  emissivities  computed  using  Ref.  17. 

Step  4:  radiation  computation.  The  UV  radiation  is  calculated  in  an  overlay  mode  using  the  flow 
solutions  obtained  at  previous  steps.  A  parallel  Monte  Carlo  radiation  code,  NEMO,  developed  at  ERC, 
that  has  2D  and  3D  modules  and  generally  allows  for  prediction  of  radiation  from  two-phase  rocket  plumes 
in  the  UV,  visible,  and  IR  ranges,  was  used.  To  calculate  the  emissivity  of  alumina  and  soot  particles,  it  is 
necessary  to  know  optical  properties  of  the  corresponding  materials.  The  real  part  of  the  complex  refraction 
coefficient  of  alumina  particles  was  calculated  using  the  expression  given  in  Ref.  18.  The  imaginary  part  was 
calculated  using  two  different  sets  of  expressions,  (i)  set  of  Ref.  16,  and  (ii)  set  of  Ref.  17.  When  calculating 
the  absorption  index  according  to  Ref.  17,  to  provide  for  a  smooth  transition  between  the  fundamental 
absorption  and  the  Urbach  edge  absorption,  a  term  bXopt  was  added  to  the  expression  for  the  fundamental 
absorption,  similar  to  Ref.  16.  For  soot,  the  complex  refraction  index  was  calculated  according  to  Ref.  19. 
Once  the  optical  properties  of  the  materials  in  the  temperature  and  wavelength  ranges  of  interest  (300- 
3,000  K  and  0.2-20  pm,  respectively)  were  calculated,  the  absorption  efficiency  Qabs  was  calculated  for  the 
given  particle  radii  using  the  BHMIE  code  first  published  in  Ref.  20.  Then,  the  particle  emissivity  e  was 
obtained  as 

^ max 

e=  I  QabsEbxdX/aT 4  (1) 

A  min 

where  Eb\  is  the  black  body  emissive  power,  a  is  the  Stefan-Boltzmann  constant,  A mjn  =  0.2  pm,  and 

A  max  ~  20  /irn. 

The  two  different  sets  of  expressions  result  in  significant  differences  in  the  absorption  index  for  the 
range  of  temperatures  important  in  this  work,  as  illustrated  in  Fig.  3  (left).  Note  that  the  newer  data 
set  is  characterized  by  much  lower  absorption  index  for  higher  wavelengths,  which  is  expected  to  strongly 
impact  the  UV  spectra.  The  difference  is  also  reflected  in  particle  emissivities  computed  with  these  sets. 
Figure  3  (right)  shows  emissivities  for  0.2  and  3.66pm  diameter  particles  as  a  function  of  surface  temperature, 
calculated  using  Ref.  16  and  Ref.  17.  As  expected,  the  emissivity  is  much  higher  for  bigger  particles. 
The  temperature  dependence  is  similar  for  both  diameters.  For  the  model  of  Ref.  16,  the  emissivity  stays 
approximately  constant  until  the  temperature  reaches  ~1750  K,  and  then  grows  rapidly.  For  the  model  of 
Ref.  17,  the  emissivity  starts  growing  after  2000  K.  In  the  2000-3000  K  temperature  range,  the  emissivity 
predicted  by  model16  is  much  higher  than  both  liquid  and  7-phase  emissivities  predicted  by  the  model.1' 

6  of  21 

American  Institute  of  Aeronautics  and  Astronautics 
Distribution  A:  Approved  for  public  release;  distribution  unlimited 


IV.  Results  and  Discussions 


A.  Gas  Properties  in  the  Plume 

The  flow  fields  that  illustrate  gas  properties  from  the  nozzle  exit  to  2.5  km  downstream  are  given  in  Fig.  4. 
The  three  overlapping  domains  here  illustrate  the  multi-zone  DSMC  approach  used  in  this  work  (each 
successive  zone  uses  the  solution  from  the  previous  zone  for  the  inflow  boundary  condition) .  The  translational 
temperature  field  clearly  shows  the  increase  in  temperature  in  the  mixing  layer  between  the  plume  and  the 
free  stream.  Note  that  the  25°  line  of  sight  of  the  far  field  photometer  intersects  with  the  elevated  temperature 
region  at  about  500  m  from  the  nozzle.  The  translational  temperature  in  that  region  is  about  1,000  K,  which 
should  not  have  any  implications  for  UV  radiation,  except  for  possible  evaporation  of  condensed  liquid 
droplets  (that  are  not  accounted  for  in  this  work).  Typical  for  plumes  expanding  in  a  background  gas,  the 
compression  region  produced  by  the  interaction  of  the  free  stream  and  the  plume  converges  to  the  nozzle 
axis  at  some  distance  from  the  nozzle.  In  this  case,  this  distance  is  over  1  km,  as  the  Mach  number  field 
clearly  shows.  The  Mach  number  at  the  nozzle  axis  reaches  its  minimum  of  about  5  at  about  1,400  m  from 
the  nozzle.  The  maximum  Mach  number  of  about  30  is  observed  in  the  expanding  plume  at  X«600  m. 

Although  the  mean  free  path  in  the  free  stream  is  large,  and  the  Knudsen  number  based  on  the  rocket 
length  of  3.1  m  is  about  1,  the  flow  in  the  mixing  layer,  which  is  hundreds  of  meters  wide  and  kilometers  long, 
is  near-continuum.  This  is  illustrated  in  Fig.  4  (bottom)  where  the  gas  mean  free  path  is  plotted.  Generally, 
this  makes  conventional  continuum  solvers  applicable  in  that  region.  However,  there  is  an  expansion  region 
in  the  core  flow  where  the  gas  mean  free  path  continues  to  increase,  reaching  its  maximum  of  about  40  m. 
The  radius  of  the  plume  core  flow  varies  from  100  to  200  m,  which  is  comparable  to  the  gas  mean  free  path 
there.  The  applicability  of  continuum  approaches  in  that  region  may  therefore  be  questionable. 

B.  Alumina  Particle  Properties  in  the  Far  Field 

Comparison  of  gas  flow  fields  obtained  for  different  alumina  particle  emissivities  and  particle  size  distributions 
have  shown  that  the  effect  of  these  properties  on  gas  in  the  far  field  of  the  plume  is  relatively  small.  The 
maximum  difference  observed  between  Cases  1  and  4  amounted  to  about  5  percent  for  gas  temperature,  and 
still  smaller  for  gas  densities.  As  expected,  the  largest  effect  of  particle  properties  is  seen  in  particle  surface 
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temperatures.  The  comparison  of  the  most  important  property  in  terms  of  UV  radiation,  namely,  the  surface 
temperature  of  3.66/im  alumina  particles,  is  presented  in  Fig.  5. 

The  upper  two  plots  illustrate  the  impact  of  the  particle  enrissivity.  The  values  of  the  enrissivity  in  Case  1 
are  up  to  an  order  of  magnitude  larger  than  those  in  Case  2.  This  results  in  much  faster  cooling  of  particles. 
For  Case  1,  the  particle  surface  temperature  decreases  to  the  critical  value  of  2,100  K  at  55  m  from  the 
nozzle  exit.  For  Case  2,  it  takes  particles  over  400  m  to  reach  this  temperature.  After  the  crystallization 
begins,  for  both  cases  the  temperature  reaches  the  alumina  melting  temperature  of  2,320  K.  This  value  is 
reached  at  about  500  m  for  Case  2  instead  of  60  m  for  Case  1.  One  may  expect  larger  effects  of  the  phase 
change  for  Case  1  than  for  Case  2,  since  particle  densities  are  very  small  at  500  m.  The  influence  of  the  size 
distribution  on  the  particle  surface  temperature  is  negligibly  small  compared  to  that  of  particles  emissivities 
for  both  large  and  small  particles  (compare  Fig.  5  middle  and  bottom).  Case  4  is  practically  identical  to 
Case  3,  and  is  not  shown  here. 

Finally,  the  density  of  sub-micron  particles  in  the  far  field  is  presented  in  Fig.  6  for  Case  2  (the  other  cases 
are  qualitatively  similar).  As  follows  from  the  near-field  results  presented  above,  there  is  a  large  number  of 
sub-micron  particles  in  the  field  of  view  of  the  25°  photometer.  This,  along  with  high  temperature  of  large 
particles  in  the  far  held,  means  that  the  contribution  of  the  far  held  sub-micron  particles  on  the  radiance 
detected  by  the  25°  instrument  may  be  important  and  needs  to  be  analyzed.  However,  the  temperature 
of  these  particles  does  not  change  much  in  the  far  held,  and  stays  at  around  1,000  K,  which  is  too  low  for 
signihcant  UV  emission. 

C.  Soot  Particle  Properties  in  the  Far  Field 

The  properties  of  soot  in  the  plume  are  now  examined.  Although  only  Case  1  is  shown  here,  the  soot 
properties  for  the  other  three  cases  are  very  similar.  The  number  density  of  100  nm  soot  particles  is  shown 
in  Fig.  7  (left).  Similar  to  sub-micron  alumina,  a  large  amount  of  soot  reaches  the  held  of  view  of  the  25° 
photometer.  Light  and  small  soot  particles  to  the  large  extent  follow  the  gas  streamlines  in  the  near  held, 
although  then  deviate  slightly  from  gas  streamlines  in  the  far  held.  The  gas  in  the  far  held  is  too  rarefied  to 
signihcantly  affect  the  soot  temperature,  although  there  is  a  visible  increase  in  temperature  at  larger  angles 
from  the  plume  axis,  where  the  soot  particles  interact  directly  with  the  free  stream.  Even  though  the  soot 
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emissivities  are  much  higher  than  those  of  same-sized  alumina  particles,  low  soot  temperatures  preclude  soot 
from  being  direct  contributor  the  the  UV  emission;  only  scattering  processes  are  important. 

The  number  density  profiles  of  soot  particles  along  the  line  of  sight  of  the  25°  photometer  are  given  in 
Fig.  8.  Very  few  particles  are  near  the  instrument.  Note  here  that  the  values  in  the  first  few  centimeters 
may  actually  be  even  lower  if  soot  particles  were  modeled  inside  the  nozzle,  and  not  just  assumed  to  follow 
the  gas  streamlines  there.  Still,  it  is  believed  to  be  a  minor  factor  in  terms  of  the  actual  contribution  of 
soot  particles  in  photon  scattering,  especially  keeping  in  mind  the  uncertainties  in  the  soot  formation  and 
distribution.  The  soot  density  reaches  its  maximum  at  about  8  m  from  the  detector,  and  then  decreases 
almost  linearly  (in  log  scale)  in  the  far  field.  There  is  some  deviation  from  linearity  for  smaller  particles, 
that  is  attributed  to  the  impact  of  the  free  stream  on  soot  particle  trajectories. 

D.  UV  Radiance  at  4°  and  25° 

The  availability  of  experimental  data5  allows  quantitative  analysis  of  the  impact  of  different  factors  mentioned 
above  on  the  radiation  intensity.  The  spectral  radiance  at  230  nm  is  presented  in  Table  2  for  four  numerical 
cases  under  consideration,  along  with  the  experimental  results.  Two  photometer  lines  of  sight  are  considered, 
4°  and  25°.  Consider  first  the  the  near  field  4°  results.  A  significant  decrease  in  particle  emissivities, 
illustrated  in  the  comparison  of  Cases  1  and  2,  results  in  an  about  20%  increase  in  radiance.  This  is  related 
primarily  to  the  fact  that  the  temperature  of  large  (micron-sized)  particles  in  the  plume  decreases  much 
slower  for  Case  2.  The  higher  temperatures  cause  more  photons  to  be  emitted,  even  though  the  emissivities 
are  lower.  It  is  interesting  to  note  that  there  is  a  further  increase  in  the  radiation  when  more  sub-micron 
particle  bins  are  included  in  the  computations  (compare  Cases  2  and  3).  The  reason  for  this  is  photon 
scattering  on  small  particles.  The  large  number  of  0.2/xrn  particles  in  Case  2  results  in  a  significant  reduction 
in  radiation  due  to  scattering.  When  the  radiation  is  computed  without  sub-micron  alumina  particle  bins, 
the  radiance  values  for  Cases  2  and  3  are  almost  the  same  (8.3TCV6  W//xm  m2  sr  for  Case  3  against  8.4-KV6 
W//zm  m2  sr  for  Case  2),  since  the  number  density  and  surface  temperature  of  large  particles  for  these  two 
cases  are  almost  identical.  The  larger  number  of  micron-sized  particles  and  small  number  of  submicron-sized 
particles  result  in  a  smaller  impact  from  scattering  and  the  highest  radiation  for  Case  4. 

All  numerical  values  for  the  4°  photometer  are  noticeably  larger  than  the  measured  radiation  of  3.6T0-6 
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W//zm  m2  sr.  This  is  primarily  attributed  to  the  impact  of  the  gas-particle  heat  transfer  model.  Previous 
analysis21  have  shown  that  the  heat  transfer  model  implemented  in  CFD++  is  reproduced  closely  by  the 
free  molecular  model  with  the  energy  accommodation  coefficients  of  about  0.4.  A  model  with  more  efficient 
heat  transfer  (larger  energy  accommodation  coefficient)  will  cool  particles  inside  the  nozzle  more  effectively, 
and  thus  the  radiative  emission  would  decrease.  Another  reason  is  the  likely  underprediction  of  the  number 
of  sub-micron  particles  in  the  numerical  simulation,  as  discussed  below. 

The  result  given  in  Table  2  show  that  the  radiation  computed  along  the  25°  line  of  sight  is  significantly 
lower  than  the  measured  value.  The  maximum  radiation,  recorded  for  Case  3,  is  almost  three  times  lower 
than  the  data.  As  all  the  radiation  comes  to  the  25°  detector  due  to  the  scattering  on  sub-micron  particles,1 
it  is  clear  that  the  model  underpredicts  the  amount  of  these  particles  in  the  plume.  The  formation  of  a 
somewhat  larger  number  of  sub-micron  alumina  particles  is  possible,  as  well  as  larger  soot  concentrations 
in  the  exhaust.  Note  also  that  the  current  work  included  only  two  out  of  four  possible  sources  of  sub¬ 
micron  particles,  aluminum  oxidation  and  soot  formation.  The  other  two,  heterogeneous-  and  homogeneous 
condensation  and  nozzle  surface  degradation,  are  not  modeled. 

An  assessment  of  condensation’s  potential  impact  on  the  results  was  conducted.  The  calculation  of  the 
water  supersaturation  ratio  has  shown  that  the  condensation  will  start  at  about  10  m  from  the  nozzle  exit 
for  25°  line  of  sight  and  20  nr  for  4°  line  of  sight,  because  of  the  significant  decrease  in  gas  temperature.  The 
analysis  of  collision  rates  in  the  plume  provides  the  upper  estimates  for  the  water  cluster  sizes  in  the  plume 
(sticking  is  assumed  to  be  unity  and  no  evaporation  is  allowed),  which  are  about  1,200  monomers  per  cluster. 
It  is  clear  that  the  size  of  such  clusters  is  too  small  to  be  a  significant  factor  in  UV  photon  scattering.  The 
heterogeneous  condensation  on  HC1  molecules  is  also  not  expected  to  be  a  major  factor  for  UV  radiation 
signatures.  The  impact  of  another  source  of  submicron  particles,  nozzle  surface  degradation,  is  expected  to 
be  well  below  one  percent  due  to  their  relatively  small  number. 

The  impact  of  soot  particles  in  the  plume  on  the  spectral  radiation  at  230nm  is  shown  in  Table  3.  The 
results  indicate  that  the  inclusion  of  0.3%  soot  results  in  a  25%  increase  in  radiation  in  the  far  field  for  Case 
2  (the  results  for  Case  3  are  similar,  whereas  for  Case  4  the  photons  reflected  on  soot  particles  constitute 
about  two  thirds  of  the  total  radiation  at  25°).  It  is  important  to  note  that  all  three  bins  of  soot  particles 
(even  the  smallest  bin  of  40  nnr)  contribute  noticeably  to  the  far  field  radiation.  The  increase  of  soot  mass 
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loading  in  the  plume  from  0.3  to  3%  results  in  a  drastic  increase  in  the  far  field  radiation,  with  the  25°  line 
of  sight  radiation  being  only  about  20%  lower  than  the  data  (although  such  a  high  mass  loading  of  soot  does 
not  seem  very  likely).  The  near  field  radiation,  that  does  not  change  much  for  the  0.3%  soot  case,  decreases 
about  2%  for  the  3%  soot  case. 

In  addition  to  the  uncertainties  related  to  the  gas-particle  heat  transfer  and  particle  size  distribution, 
that  both  influence  the  agreement  between  the  experimental  and  numerical  data,  there  is  also  an  uncertainty 
related  to  the  accuracy  of  the  line  of  sight  value  determination.  The  actual  value  may  somewhat  differ  from 
the  reported  4°  and  25°,  and  a  small  difference  in  the  line  of  sight  results  in  a  noticeable  difference  in  the 
radiation  prediction,  as  indicated  in  Table  4.  Relatively  large  numbers  of  submicron  particles  in  the  far 
field  raise  the  question  of  convergence  of  the  obtained  results  in  terms  of  the  length  of  the  computational 
domain.  In  order  to  prove  the  convergence,  the  radiation  calculations  have  been  conducted  with  varying 
domain  length  (in  all  runs  only  the  downstream  boundary  was  moved).  The  results  for  the  25°  line  of  sight 
are  presented  in  Fig.  9.  The  main  conclusion  here  is  that  the  radiation  shows  minimum  change  after  the 
first  few  hundred  meters  downstream  from  the  nozzle  axis.  Note  that  the  near  field  radiance  converges  even 
faster.  That  means  that  the  DSMC  modeling  of  the  plume  near  and  mid-fields  only,  without  the  far  field, 
would  provide  adequate  accuracy  for  radiation  predictions. 

Previous  results1  have  shown  a  qualitative  difference  in  the  shape  of  the  UV  spectra  from  the  published 
data.  While  the  agreement  was  reasonable  (within  a  factor  of  two)  for  200-300  nm  range,  at  400  mn  the 
numerical  results  overpredicted  by  a  factor  of  20.  This  was  attributed  mostly  to  the  influence  of  alumina 
particle  emissivities  (the  absorption  index  of  Ref.  16  was  used  in  Ref.  1).  Similar  to  the  previous  results,  the 
computed  spectral  shape  is  much  steeper  than  the  experimental  one  when  the  emissivities  based  on  Ref.  16 
are  used,  as  illustrated  in  Fig.  10,  Case  1.  The  emissivities  computed  from  Ref.  17,  however,  allow  one  to 
capture  the  shape  of  the  experimental  spectra  much  closer  (Cases  2  through  4). 

V.  Summary  and  Conclusions 

A  computational  study  of  a  Star  27  solid  propellant  rocket  plume  at  an  altitude  of  118  km  has  been 
conducted  using  a  combined  multi-step  continuum/kinetic  approach.  The  solution  of  the  Navier-Stokes 
equations  was  obtained  inside  the  nozzle  and  in  the  plume  in  the  vicinity  of  the  nozzle  exit.  The  DSMC 
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method  was  used  to  simulate  most  of  the  plume,  as  well  as  the  plume-free  stream  interaction.  A  Monte 
Carlo  based  radiation  solver  was  used  in  an  overlay  mode  to  calculate  UV  radiation  in  the  plume.  All  gas 
dynamic  flow  solvers  used  in  this  work  included  two-way  coupling  between  the  gas  flow  and  particles,  with 
two  types  of  particles  being  considered,  alumina  and  soot. 

The  impact  of  alumina  particle  emissivities  on  particle  properties  and  UV  radiation  in  the  near  and  far 
field  has  been  analyzed.  It  has  been  shown  that  the  use  of  emissivities  based  on  Ref.  17  results  in  much  lower 
particle  temperatures  in  the  plume  and  significant  delay  of  the  liquid-to-7  phase  change,  as  compared  to 
the  case  when  emissivities  based  on  Ref.  16  are  used.  The  radiation  was  found  to  be  tens  of  percent  higher 
both  in  the  near  and  far  field  when  smaller  values  of  emissivities  are  used  due  to  slower  particle  cooling. 
The  effect  of  particle  emissivities  on  the  gas  is  relatively  small.  Comparison  with  available  experimental 
data  have  shown  that  the  model  overpredicts  the  experimental  data  in  the  near  field  (4°  line  of  sight)  by  on 
average  a  factor  of  two,  which  is  thought  to  be  the  effect  of  the  gas-particle  heat  transfer  model  inside  the 
nozzle. 

The  numerical  results  for  the  far  field  (25°  line  of  sight)  UV  radiation  underpredict  the  measurements, 
with  the  difference  ranging  from  almost  a  factor  of  ten  to  about  25%,  depending  on  the  sub-micron  alumina 
particle  size  distribution  and  soot  mass  loading.  Both  these  factors  are  not  well  known,  and  were  therefore 
parameterized.  Hot  micron-sized  particle  emission  in  the  plume  core  flow  scattered  on  sub-micron  particles 
in  the  side  flow  is  the  dominant  source  of  the  far  field  UV  radiation,  and  both  the  alumina  particles  and 
soot  were  found  to  significantly  contribute  to  the  scattering  process.  The  effect  of  condensation  and  nozzle 
surface  erosion  on  the  UV  photon  scattering  was  estimated  to  be  relatively  small.  Both  far  and  near  field 
radiation  prediction  were  found  to  be  sensitive  to  the  accuracy  of  the  line  of  sight  determination. 
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Table  1.  Mole  Fractions  (M.F.)  of  major  plume  species. 


Species 

CO 

C02 

Cl 

H 

h2 

h2o 

HC1 

OH 

n2 

M.F. 

0.2225 

0.0291 

8.67xl0~3 

0.0190 

0.2734 

0.1767 

0.1590 

5.38xl0-4 

0.1108 

Table  2.  Radiation  at  230  nm  for  different  models,  W //mi  m2  sr.  0.3%  soot  is  included. 


Case 

4° 

25° 

Case  1 

4.85-10"6 

3.55-10-9 

Case  2 

5.75-10-6 

6.20-10"9 

Case  3 

6.99-10-6 

7.70-10-9 

Case  4 

9.45-10"6 

2.28-10-9 

Experiment 

3.6-10-6 

22-10"9 
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Table  3.  Effect  of  soot  particles  on  radiation  at  230  nm,  Case  2,  W//zm  m2  sr. 


Line  of  sight 

Soot  content 

Radiance 

25° 

no  soot 

4.83-10"9 

25° 

40  nm  bin  only 

5.21-10-9 

25° 

40  nm  and  100  nm  bins 

5.92-10-9 

25° 

all  soot  bins 

6.20-10-9 

25° 

3%  soot 

17.7-10-9 

4° 

no  soot 

5.75-10-6 

4° 

3%  soot 

5.62-10-6 

Table  4.  Effect  of  the  line  of  sight  value  on  radiation  at  230  nm,  Case  2,  W//im  m2  sr. 


3.0 

7.43-10"6 

24.0 

9.18-10-9 

3.5 

6.57-10-6 

24.5 

7.60-10-9 

3.9 

5.89-10-6 

24.9 

6.47-10-9 

4.0 

5.75-10-6 

25.0 

6.20-10-9 

4.1 

5.58-10-6 

25.1 

5.95-10-9 

4.5 

5.07-10-6 

25.5 

5.06-10-9 

5.0 

4.60-10"6 

26.0 

4.09-10-9 

16  of  21 


American  Institute  of  Aeronautics  and  Astronautics 
Distribution  A:  Approved  for  public  release;  distribution  unlimited 


icr  io0  To1 

Diameter,  |im 


Figure  2.  Three  particle  size  distributions. 
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Figure  3.  Alumina  absorption  index  calculated  with  two  different  sets  of  expressions,  Ref.  16  and  Ref.  IT  (7 
phase),  for  a  surface  temperature  of  2,200  K  (left).  Emissivities  for  two  particle  diameters  computed  with 
these  two  sets  (right). 


X,  m 


2400 

1821 

1382 

1048 

795 


30.00 

19.02 

12.06 


40.00 

20.56 

10.56 


Figure  4.  Plume  translational  temperature  in  Kelvin  (top),  Mach  number  (middle),  and  gas  mean  free  path 


in  meters. 
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Figure  5. 


Surface  temperature  field  of  3.66/im  particles  for  Case  1,  top;  Case  2,  middle;  Case  3,  bottom. 


X,  m 


1.0E+11 
4.4E+10 
8.5E+09 
3.2E+08 
1 .4E+08 
6.1E+07 
2.7E+07 
1 .2E+07 
5.2E+06 
2.3E+06 
1.0E+06 


Figure  6.  Number  density  field  in  kg/m3  of  0.2/im  particles. 


X,  m 


1.4E+11 

1.8E+10 
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1.6E+07 

6.0E+06 

2.2E+06 
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X,  m 


5.00E+02 

4.80E+02 
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Figure  7. 


Number  density  in  particle/m3  (left)  and  surface  temperature 


in  Kelvin  of  0.1/irn  diameter  soot 


particles. 
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Figure  8.  Soot  number  density  along  the  25°  line  of  sight. 


Figure  9.  The  impact  of  the  far  field  on  the  total  radiance  at  230  nm,  Case  2. 
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Figure  10.  Comparison  of  spectral  radiance  for  different  cases  with  experimental  data. 
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